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We investigate interplay between magnetic fluctuations and superconductivity in the effective 
five-band Hubbard model for iron-oxypnictide superconductors on the basis of the fluctuation- 
exchange approximation. As for the normal-state properties, we find the pseudogap behavior in 
the NMR relaxation rate and the spectral weight in the electron-doped region, while we cannot 
find such behavior in the hole-doped region. The pseudogap behavior originates from the band 
structure effect, that is, existence of high density of states just below the Fermi level. 

Solving the superconducting Eliashberg equation, we find that the most probable candi- 
date for the pairing symmetry is the sign-changed s-wave spin-singlet state. For small Hund's 
coupling J, the eigenvalue is not so sensitive to carrier doping, and seems to be irrelevant 
with antiferromagnetic (AF) spin fluctuation. We suggest that correlation between spin and 
spin-quadrupole is important as the pairing mechanism as well as the AF fluctuation. 
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The recent discovery of the Fe-based layered supercon- 
ductor LaFeAsOi-^F^ 1 with the transition temperature 
T c = 26K has promoted highly intensive research activi- 
ties. In the new family of FeAs superconductors, T c has 
been elevated to over 50K by substituting rare earth for 
La, 2 which is the highest after the high-T c cuprates. 

Many experimental and theoretical efforts have re- 
vealed the following common features. The undoped sys- 
tem undergoes the stripe-type antiferromagnetic (AF) 
transition with the structural phase transition. 3 ' 4 With 
electron doping, the system shows the superconducting 
transition. 1 The hole-doped system 5 7 also have the same 
order of T c , which is insensitive to carrier doping and the 
AF fluctuation. 1,4,5 In the band calculation, 8 the disper- 
sion near the Fermi level is mainly composed of Fe 3d 
orbitals, which are occupied by six electrons in the un- 
doped system. The Fermi surfaces (FSs) are composed of 
a 3-dimensional one around Z, two hole cylinders around 
T and two electron cylinders around M. The dispersion 
relation and the FSs have been roughly verified by the an- 
gle resolved photoemission spectroscopy (ARPES). 9 As 
for the superconducting gap symmetry, the NMR Knight 
shift indicates the spin-singlet state. 10 14 T c seems to 
be robust for impurities. 14 The penetration depth, 15 the 
specific heat 16 and ARPES 17, 18 indicate almost isotropic 
two gap system, while the NMR relaxation rate 1/Xi 
shows no coherence peak and the T 3 law, 4, 10 13, 19 which 
suggest existence of line nodes. At present, these facts 
are controversial. 20 Conventional theory by the electron- 
phonon interaction cannot explain high T c . 21 We can 
expect that unconventional superconductivity originates 
from the electron correlation. 

Among many theoretical research studies, the most 
promising pairing state is the sign-changed s-wave spin- 
singlet state (hereafter, ,s ± -wave state). 22,23 Kuroki et 
al. 23 have indicated that the dispersion near the Fermi 
level can be well explained by simple two-dimensional 



square lattice of an Fe atom, although two Fe atoms 
are contained in the actual unit cell. Using the random 
phase approximation (RPA), they investigated the effec- 
tive five-band Hubbard model with the Coulomb inter- 
action on an Fe site; the intra-orbital Coulomb [/, the 
inter-orbital Coulomb U', the Hund's coupling J and the 
pair-hopping term J' . Recently, by the third-order per- 
turbation theory, 24 and RPA in 16-band d-p model, 25 the 
possibility of such s^-wave state has been verified within 
weak coupling approach. In this paper, we investigate in- 
terplay between magnetic fluctuations and superconduc- 
tivity in the above effective five-band Hubbard model on 
the basis of the strong coupling approach. 

First, we obtain a similar band structure by the tight- 
binding fitting using hopping parameters in Table I of 
Ref.23. The band dispersion and the FSs in the unfolded 
Brillouin zone (BZ) are illustrated in Fig.l. V (n, 7r) be- 
comes equivalent to T (0,0) in the folded original BZ. 
Electron pocket around M and hole pocket around T 
(T f ) reproduce the FSs in the original band structure, 
although the FS around T f is somehow larger than that 
in the band calculation by Kuroki et al. The yz, zx and 
x 2 — y 2 orbitals have large weight near the Fermi level. 

We here examine the correlation effect for this band 
structure within the self-consistent second-order pertur- 
bation (SC-SOPT), and the fluctuation-exchange ap- 
proximation (FLEX). The FLEX calculation in the 
multiband system follows Ref.26. In the present five- 
band model, the normal (anomalous) Green's function 
Gim(k) {J-£m(k)) and self-energy Ui m (k) can be treated 
as a 5 x 5 matrix in orbital indices, while the spin and 
the charge (orbital) susceptibilities (x^/ mm /(#)) as a 
5 2 x 5 2 matrix. With x° simple bubble diagram of (5, 
X s,c — X° =t x°V S}C x s,c , where V s,c denote bare vertices 
for the spin and charge sectors. 27 We set X s = X c = X° m 
SC-SOPT, and E = in RPA. Actual numerical calcula- 
tions have been carried out mainly with 64 x 64 /c-point 
meshes and 256 Matsubara frequencies. Some calcula- 
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Fig. 1. The band dispersion (a) and the FS (b) at undoped n = 
6.00 in the unfolded BZ. In (a), red and green bars represent 
the weight of z 2 and x 2 — y 2 orbitals, respectively. In (b), red 
line around T (r x ) and green line around M represent hole and 
electron FSs, respectively. Gray line denotes the original BZ. 



tions have been checked with 64 x 64 x 512. We hereafter 
set U' = U — 2 J and J' = J in the Coulomb interaction, 
and a unit of energy to be electron volt (eV). 




Fig. 2. The band dispersion (a) and the FS (b) at n = 6.00 for 
(U,J,T) = (1.32,0.22,0.01) in SC-SOPT. 



In Fig. 2, we illustrate the band dispersion and the FSs 
included only the self-energy shift Z^ m (fe,i7rT) in SC- 
SOPT. Compared with Fig.l, z 2 orbital near r' shifts 
upward, and x 2 — y 2 orbital downward. With large J, 
these shifts increase, and a new FS with z 2 orbital ap- 
pears. Also in simple SOPT and FLEX, the tendency 
is the same. At a glance, this seems to explain a large 
outer FS observed in ARPES. 9 ' 17 ' 18 In ARPES, how- 
ever, a flat band, which is assigned to a z 2 orbital, is 
observed below the Fermi level. In FLEX, furthermore, 
with appearance of the new FS, the dominant magnetic 
fluctuation changes from (tt,0) to (7r,7r). Thus, the self- 
energy shift by the correlation effect cannot be so large 
to reproduce experimental results. This indicates that 
the Hund's coupling J should not be so large. 28 The 
large outer FS observed in ARPES is probably composed 
of the x 2 — y 2 orbital. In addition, existence of the FS 
works in favor of superconductivity. Thus, it will be nat- 
ural that the outer FS composed of x 2 — y 2 orbital ex- 
ists above z 2 orbital around V. The self-energy shift in 
question works to equalize electron number included in 
each orbital. Although this fact seems to be physically 
correct as the correlation effect, we have to investigate 
magnetic and superconducting properties without chang- 
ing the character of the FSs from the band calculation. 
This trouble is because the band calculation already in- 
cludes partially exchange-correlation effect. We should 
not overcount the static correlation effect at uo = in 



the normal self-energy. The final result which includes 
the correlation effect should reproduce the band charac- 
ter obtained by band calculations, which coincides with 
experimental results. At present, we hereafter shift site 
energies of z 2 and x 2 — y 2 orbitals by +0.1 and +0.3, 
respectively. 29 This tentative method allows us to inves- 
tigate superconductivity without changing the character 
of FSs and the dominant magnetic fluctuation. However, 
we need further investigation to proceed to more quan- 
titative explanation. 
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Fig. 3. (a) Partial DOS and (b) temperature dependence of total 
DOS at electron-doped n = 6.12 for (U, J, T) = (1.44, 0.24, 0.01). 



Now, let us show numerical results in FLEX with site 
energies shifted. In Fig. 3, we display the partial DOS 
and its temperature dependence. The yz, zx and x 2 — y 2 
orbitals hold large weight at the Fermi level, and high 
DOS just below the Fermi level. Corresponding to this 
fact, we find the pseudogap behavior in Fig. 3(b). On the 
contrary, we cannot find such behavior in the hole-doped 
region, although we do not demonstrate it. These be- 
haviors are consistent with the photoemission spectra, 30 
qualitatively. Thus, this pseudogap can be naturally ex- 
plained from the band structure. This is the same reason 
as that proposed by Yada et al. for Nao.35Co02. 31 

Next, let us evaluate the NMR relaxation rate 1/Ti, 
given by 1/T\T oc Imxf oc (^)/o;| w . Using the Pade 
approximation, we carry out numerically analytic con- 
tinuation for the local spin susceptibility, xf oc = 
X^, m Xi\ mm W) obtained in FLEX. As illustrated in 
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Fig. 4. NMR 1/TiT at intervals of An = 0.04 for (U, J) 
(1.44,0.24). The pseudogap appears in the electron-doped. 



Fig. 4, with decreasing temperatures, \jT\T is enhanced 
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in the undoped and hole-doped systems. On the other 
hands, in the electron-doped system, it shows the pseu- 
dogap behavior. It becomes remarkable with electron- 
doping. This qualitatively explains the experimental re- 
sult in NMR. 4 ' 11 13 ' 19 ' 32 Such pseudogap behavior ex- 
ists even without electron correlation in the present band 
structure. The electron correlation only reduces the en- 
ergy and temperature scale to actual values. 





i 




0.8 




0.6 








0.4 




0.2 








°0 



n=5.92 



n=6.12 



(a) 


1 

T=0.080 

T=0.020 






• T=0.008 























■ 

T=0.080 

v T=0.020 




/ * 




- T=0.008 



























3 



Fig. 5. Imxf oc (<*0 at hole-doped n = 5.92 (a) and electron-doped 
n = 6.12 (b) for (J7, J) = (1.44, 0.24). 



In Fig. 5, we show temperature dependence of imagi- 
nary part of the local spin susceptibility, Imxf oc (o;). The 
slope at oo = corresponds to 1/T\T. We find that differ- 
ences between the electron- and the hole-doped systems 
are restricted to the low frequency region, and the peak 
position at oo ~ 1 is almost unchanged, although the 
spectral weight is broad in the hole-doped. This fact may 
be related to weak sensitivity of T c for carrier doping. 

Next, we solve the superconducting Eliashberg equa- 
tion for the spin-singlet pairing within FLEX. The most 
probable candidate is the s^-wave state. The second 
largest eigenvalue is given by d x i _ y 2 -wave. In Fig. 6, we 




of the gap function around T f is about half. 34 This is 
consistent with magnitude of the gap on the outer FS 
obtained in the ARPES analysis. 17 ' 18 
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Fig. 7. Phase diagram for J = 0.22 (a) and 0.24 (b) at U 
1.44. 35 A denotes the eigenvalue for the s^-wave. 



In Fig. 7, we illustrate contour map of eigenvalues 
for the ,s ± -wave state and the static spin susceptibility 
X s (Q,0) at Q = (tt,0) for J = 0.22 and 0.24. With 
carrier doping, the AF spin fluctuation decreases, while 
eigenvalues for J = 0.22 are not so sensitive. Eigenvalues 
for J = 0.24 have stronger correlation with the AF spin 
fluctuation. 35 The phase diagram for small J seems to 
be consistent with experiments, 1,5 although we need to 
evaluate T c itself actually. This point is related to the 
pairing mechanism discussed below. 

Since the early stage, the pairing mechanism has prob- 
ably been thought to be AF spin fluctuation, 22,23 which 
is related to the AF state in the undoped system. 3 Surely, 
the AF spin fluctuation is important. However, this can- 
not explain insensitivity of T c for carrier doping 1,5 and 
weak correlation with the AF fluctuation. 4, 13 We need 
to make clear the pairing mechanism. 

We hereafter study the mechanism within RPA, since 
the normal self-energy is unnecessary for this analysis. 
First, let us examine spin and charge fluctuations, since 
the pairing interaction can be divided into the spin sec- 
tor, 3V s x s V s /2, and the charge sector, V c x c V c /2. In 
addition, we can carry out analysis like perturbation the- 
ory without vertex corrections by truncating x s,c within 
finite-order terms. In Table I, we indicate the result for 
the s^-wave in the linearized gap equation. Within the 



Fig. 6. Sum of the third- and fourth-band anomalous Green's 
functions at n = 6.00 for (U,J,T) = (1.44, 0.24, 0.008). 33 ' 34 It 
has about +1 on the outer FS around T, — 1 around M and +0.5 
around F' . The second-band possesses about +1 on the inner FS 
around T, and the first- and the fifth-band contributions are very 
small, although we does not display them. 

display the sum of the third- and fourth-band anomalous 
Green's functions, 33 which indicates no nodal structure 
on the FSs more clearly than the gap function itself. It 
has large weight on the outer FS around T and the FSs 
around T f and M. The sign around T (V) is opposite to 
that around M. In addition, we notice that magnitude 



Table I. Eigenvalues in the truncated perturbation at n = 6.1 
for (U,J : T) = (1.00,0.10,0.008). V s , V c and Total represent, 
respectively, eigenvalues with only spin sector, only charge sector, 
and the total contribution. Order 00 corresponds to RPA. 
Order V~ s V~ c Total 
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third- and fifth-order perturbation, the spin sector and 
the charge sector work cooperatively, and the latter has 
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dominant contribution. However, in the second- and the 
fourth-order, the charge sector works negative. Thus, 
the contribution of charge sector oscillates in each order 
term. In such a case, by summing up to infinite order 
terms like RPA, we can expect to obtain the correct re- 
sult. Consequently, the charge sector does not contribute 
to the pairing interaction. The dominant contribution 
comes from the spin sector. 

Next, let us investigate individual matrix elements of 
X s sandwiched between V s . Even if components includ- 
ing xy and z 2 orbit als are set to zero, the eigenvalue is 
almost unchanged. On the other hand, the interaction 
among yz, zx, x 2 — y 2 orbitals is especially important. 
These three components are entangled completely. 
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Fig. 8. Correlation between the eigenvalue A and some dominant 
fluctuations for (a) (U, J) = (1.0,0.2) and (b) (1.2,0.0) within 
RPA. 36 A S pi n represents its weight which can be explained by 
only the AF spin fluctuation. x sp i n = X s (tt, 97r/64), and Xtotal is 
the sum of Xspin and x sQ • Xquad represents x QQ - 3T 



In Fig. 8, we show correlation between the eigenvalue 
A and dominant AF fluctuations for J = [7/5 and 0. 
A sp in represents the contribution from only spin-spin cor- 
relation xu.mm- It can be obtained from the eigenvalue 
when any other components of yf in the spin sector are 
set to zero. For large J = U/b in Fig. 8(a), 60 ~ 70% 
of eigenvalue can be explained by only spin-spin corre- 
lation, and the pairing interaction is dominated by the 
AF spin fluctuation. In this case, the behavior of A well 
correlates with the AF spin susceptibility X s (Q — Q?0)- 
This is inconsistent with experimental results. 1,4,5, 13 On 
the contrary, for small J — in Fig. 8(b), the spin-spin 
correlation leads to nothing but ~ 30% of eigenvalue, 
and the main part of the pairing interaction comes from 
the other fluctuation. Correlation between A and the 
AF spin fluctuation is weak. Instead, the fluctuation like 
Xu £m+Xfe ml becomes important. This is the correlation 
X s ® between spin (Su) and spin-quadrupole (Qe m )- 37 In- 
dependent of J, the behavior of A well correlates with 
Xtotal = Xspin + X s ® • Since A is a measure of the pairing 
interaction, this indicates that these fluctuations work 
cooperatively as the pairing interaction. Furthermore, it 
is interesting that the spin-quadrupole fluctuation x®® 
is enhanced for J = 0, although what is more important 
as the pairing interaction is X s ® - 

In summary, we investigated the five-band Hubbard 
model for the iron-oxypnictide superconductors mainly 



by using FLEX. We found the pseudogap behavior in 
the electron-doped region. The most probable supercon- 
ducting pairing state is the sign-changed s- wave. It is 
proper that the Hund's coupling J is not so large. In this 
case, correlation between T c and the AF spin fluctuation 
is weak, and then T c will be insensitive to carrier doping, 
rather correlates with the AF spin-quadrupole fluctua- 
tion x®® • The dominant contribution to the pairing in- 
teraction comes from correlation between spin and spin- 
quadrupole x s ® as the alternative to the AF spin fluctua- 
tion. This novel pairing mechanism naturally comes from 
the fact that the system has many degrees of freedom in 
multiorbital system, and three bands which trigger un- 
conventional superconductivity are strongly entangled. 
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With the unitary matrix w^ n , which diagonalizes the unper- 
turbed Hamiltonian in the orbital representation, we can ob- 



tain the band-diagonal JF nn (/c) = 



1 u mn^~ imi}^) • Al- 



though we show the sum of the third- and fourth-band for 
want of space, they have large weight on the different FSs. 
To examine the case with the large outer FS, we here show the 
result at n = 6.00 not electron-doped region. The magnitude 
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of gap around r' changes dependent on carrier doping. It is 
~ 0.5 in the hole-doped and ~ 1 in the electron-doped. 

35) Increase of eigenvalues around hole-doped n ~ 5.80 comes 
from change of the Fermi surface due to the self-energy shift. 
It probably does not reflect the real fact in this system. 

36) In both cases, the s^-wave state has the largest eigenvalue. 
At a glance, this seems to be inconsistent with the result in 
Ref.25. This is probably because of different dispersion rela- 
tions, especially, size of the FS around V . 

37) For yz (2) and zx (3) orbitals, for instance, spin-quadrupole 



Qf 3 is defined as J2ka ^( c k2( j c k+q3cr+cl 3cr c k+q 2a), and then, 

X2? = Q| 3 Hz)>>, and x? 3 Q = «0i 3 (9),0! 3 (-9)»- 

In Fig.8, for simplicity, we use x sQ = x!2,23+x! 2 ,32+x! 3 ,22 + 
X32 22 a ^ Q — 97r/64), included in X23 + X23 > as ^ ne dom- 
inant fluctuation, and Xquad — 

(7r,97r/64). For small J, 
an d with £, m G yz, zx, or x 2 — y 2 become large 

at q ~ Q, and they show similar temperature dependence. 
Although orbital fluctuations like Xg^ ^ m — xJ^J m ^ are a l so 
enhanced in this case, these are more suppressed by introduc- 
ing finite J as compared with x s ® an d x®® • 



